Pinwheel stabilization by ocular dominance segregation 
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^i We present an analytical approach for studying the coupled development of ocular dominance 

(OD) and orientation preference (OP) columns. Using this approach we demonstrate that OD 
^^ ■ segregation can induce the stabilization and even the production of pinwheels by their crystalliza- 

tion in two types of periodic lattices. Pinwheel crystallization depends on the overall dominance 



O 



of one eye over the other, a condition that is fulfilled during early cortical development. Increas- 
ing the strength of inter-map coupling induces a transition from pinwheel-free stripe solutions to 
F^ . intermediate and high pinwheel density states. 
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In the primary visual cortex information is processed in a two dimensional array of mod- 
ules called orientation preference (OP) columns [1]. In many species columnar patterns con- 
tain pinwheel centers, singular points around which columns activated by different stimulus 
orientations are radially arranged like the spokes of a wheel [2] . Recent research applying 
in vivo 2-photon imaging to pinwheel centers revealed that their radial organization is laid 
down with single cell precision [3]. How pinwheels are formed during visual development 
remains unresolved. In theoretical models pinwheels are generated by spontaneous symme- 
try breaking but are often dynamically unstable [4]. Recent theoretical studies, treating 
the system of orientation columns essentially as an isolated system, have examined if pin- 
wheels may be stabilized by long-range intracortical interactions [5[l, by a coupling to the 
large scale map of visual space p, or by wiring length constraints [r\. In the visual cortex, 
however, orientation columns are presumably interacting with e.g. ocular dominance (OD) 
domains, spatial frequency, and direction preference columns, see however J8[]. For instance 
OD borders intersect many of these preferentially at right angles [9]. It may thus be inad- 
equate to theoretically study the layout of orientation columns neglecting their relation to 
other columnar systems. Recently this perspective has received experimental support by a 
study reporting that orientation columns are organized more smoothly when the system of 
OD columns is removed [10]. Indeed simulations suggest that OD segregation impedes the 



process of pinwheel annihilation 



Q, 



III]. So far, however, there has been no analytic demon- 
stration that an intrinsically unstable system of orientation pinwheels can be stabilized by 
interactions with other maps. 

Here we present a dynamical systems approach for analyzing the interactions of OP and OD 
maps. We design a dynamical model for the coordinated development of OP and OD maps 
in which pinwheels become unstable in the weak coupling limit. The inter-map coupling is 
specified according to experimentally observed geometric relationships between OD and OP 
maps. Because the contralateral eye dominates during the initial formation of OD columns 



12[ | , we systematically study the impact of overall dominance by one eye on the dynamics 
of pinwheels. Using weakly nonlinear analysis we derive amplitude equations describing the 
existence and stability of pinwheel free and pinwheel rich OP maps in the coupled system. 
We identify two types of pinwheel rich solutions differing in their pinwheel density and cal- 
culate the stability and phase diagram of these solutions as a function of inter-map coupling 
and contralateral eye dominance. We find that pinwheel crystals are stable above a critical 



degree of contralateral eye dominance that induces a patchy pattern of OD domains. In- 
creasing the strength of inter-map coupling induces a transition from pinwheel free solutions 
to low and high pinwheel density patterns. In the latter regime OD segregation even induces 
the formation of additional pinwheels. 

The spatial structure of an OP map is conveniently represented by a complex field z(x) 
where x denotes the 2D position of neurons in the visual cortex, the modulus |^(x)| is a 
measure of their selectivity and ^(x) = | argz is their preferred orientation [4]. In this rep- 
resentation pinwheel centers are the zeros of the field z{-x). Ocular dominance is described 
by a real field o(x) where negative and positive values indicate ipsi- and contralateral eye 
dominance, respectively. Because OD and OP maps are not independent of each other we 
consider models containing coupling terms between both fields 

a,z(x,t) = F[z(x,t),o(x,t)] 

a,o(x,t) = G[z(x,t),o(x,t)], (1) 

where F[z, o] and G[z, ol are nonlinear operators. Various biologically detailed models have 
been cast in this form [J, ll3[]. Because cortical maps arise from a cellular instability with a 
typical wavelength A, the mathematically simplest models for the spontaneous generation 
of these patterns are of Swift-Hohenberg type [14]. We therefore choose F and G to be of 

this type and couple the fields through an energy density T 

ST' 

dt ^(x, t) = Lz ^(x, t) — \z^z — e — 

dt o(x, t) = Lo o(x, t) - o(x, tf + j-e — . (2) 

do 

Here L^o,z} = ^{o,z} ~ {^l + A) ^ ^ is an OD bias leading to an overrepresentation of the 

contralateral eye for 7 > 0, and e is the coupling strength. In this model pinwheels are 

unstable in the weak coupling limit leading to systems of stripes for e = 0, mimicking the 

behavior of competitive Hebbian models for OD or OP maps in this situation [^ . The form 

of T is found from the experimental observation that iso-orientation lines tend to intersect 

the OD borders perpendicularly [9] . T can thus be expected to contain terms of the form 

|VoV^P^. Decomposing the complex field z(x) into the selectivity \z\ and the preferred 

orientation 9 finds 

T = |VzVo|'"^ = |z|'"^(|VoVln|z||' + 4|VoV^|')'^ 

= (iVoVRezl^ + lVoVImzl^)"^. (3) 
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FIG. 1: Pinwheel annihilation and preservation in simulations of Eqn. ([2]) for different 
strengths of inter-map coupling and OD bias, Tq = 0.2, r^ = 0.02. Color coded OP map, zero 
contours of OD map superimposed, (a, b) 7 = 0,e = 0, 2000 (c, d) 7 = 0.15,e = 2000. Upper 
(lower) row: t = (lO^/r^). Initial conditions identical in (a)-(c). 



Where orientation selectivity is locally homogeneous, i.e. Vln \z\ ?^ 0, T is minimized if the 
direction of the iso-orientation lines (V^) is perpendicular to the OD borders. At pinwheel 
centers the zero contours of Re z and Im z cross and VRe z and Vim z are not parallel, T can 
be minimized only if | Vo| is small at the pinwheel centers, i.e. near extrema or saddle-points 
of o(x). In the following we analyze the case m = 2. As we will see below, this choice allows 
for a limit in which map interactions become unidirectional. 

We observe that for substantial contralateral bias and above a critical coupling e pinwheels 
are preserved or are even generated after symmetry breaking. Numerical simulations of the 
dynamics Eqn. ([2j) are shown in Fig. [H Without a contralateral bias the attractors are 
pinwheel-free stripe solutions irrespective of the strength of the inter-map coupling. 
To reveal the exact conditions for the preservation of pinwheels by inter-map coupling 
we used weakly nonlinear analysis to study the nature and stability of different types 
of solutions. To this end we first studied how the emerging OD map depended on the 
overall eye dominance. Shifting the OD field by a constant o(x, t) = o(x, t) + 5, the 
dynamics Eqn. ([2j) is mapped to 9t5(x, t) = Lo + 70^ — 5^ with L = f^ — (fc^ + A) , 
fo = To — 35^, and 7 = —36 where 6 is the real solution of —6^ + (r^ — fcc) 5 + 7 = 0; 
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FIG. 2: (a) Phase diagram of the OD model Eqn. ([2|). Dashed hnes: stabihty border of 
hexagon solutions, solid line: stability border of stripe solution, gray regions: stability region 
of homogeneous solution (b) Contralateral eye dominated neurons for the three stationary 
solutions. Circles: numerically obtained values, solid lines: Cst and C/^ex- 



an equation that has been extensively studied in pattern formation literature [15|]. It has 
three types of stationary solutions: (1) A homogeneous solution with spatially constant 



eye dominance oj^ = 5, (2) OD stripes 05t(x) = 2BstC08 {x + i/j) + 5, with Bgt = ^r/3, 
and (3) hexagonal arrays of ipsilateral eye dominance blobs in a sea of contralateral eye 
dominance Ohex{^) = lShexY.j=i^''^'^'^''' + ^•^- + ^^ ^i^h ^-fcj = 0, Y^j^j = ^^ ^^d 



^hex = —7/15 + a/(7/15) +f/15. The fractions of contralateral eye dominated terri- 
tory Cst and Chex increase with 7 as cos (CgtTr) = —S/{2Bst) and (1 — (7hex)A/3 27r ^ 
3 arccos ( ^ f — 1 + w 3 + ^^ ) ) for (2) or (3), (Fig. [Sb). The phase diagram of this model 
is depicted in Fig. [2](a). It shows the stability borders for the three solutions. Without a 
bias term the OD map is either constant, for r^ < 0, or has a stripe layout, for r^ > 0. 
For positive r^ and increasing bias term there are two transition regions. First a transition 
region from stripes to hexagons (between 7* and 72) and second a transition region from 
hexagons to the homogeneous solution (between 73 and 74). 

Close to instability, stationary solutions to the full dynamics Eqn. ([2j) can be calculated 
analytically by weakly nonlinear analysis |16|. The Fourier components of the emerging 



pattern are located on the critical circle kj = (cos jtt/S, sinjn/S) kc so that 

3 
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o(x, t) = Y. {B,{t)e^^^^ + B,{t)e-^^^') , (4) 

J 

with the complex amplitudes Aj = AjC^^^ ^ Bj = BjC^'^k Although the coupling terms 
enter at seventh order in the amplitude expansion they can be written as an effective cubic 
interaction term. Because Ai oc ^/rl and Bi oc y^, the coupling onto the OD dynamics 
becomes small for r^ <C r^, since terms like 6|A|^|5pS oc r^^rj are negligible compared to 
terms like \B\^B (x rj • In this limit, the backreaction of the OP map onto the OD map is 
thus negligible. Using uniform modes Bi = i5, the amplitude equations for the OP map are 
given by 

6 3 

Of Ai = Tz Ai — y ^ Qij \Aj\ Ai — 2 y ^ AjAj-Ai- 

6 

-ei3^5^/i,,,,A,A^A,, (5) 

with Aj- = Aj^si 9ii = ^^ 9ij = 2 and hijik an effective self-interaction tensor. The dynamics 
of the modes Ai- is given by interchanging Ai and Ai- . A solution of hexagonal symme- 
try (symmetric under rotation by 120 degree) to Eqn. ([51) is given by the uniform solution 
Aj = Aj- = A, (t)j = ^j + {j - l)27r/3 + rf5^-2, and (f)^- = -^j + {j - l)2iT/?> + d {5j^i + 5^-3), 
where we choose t/^i = 7/^3 = 0, ^02 = tt and the constant rf ?^ 1.176 is the solution of a tran- 
scendental equation. For negligible backreaction B = Bhex and A? ^ r^/ (9 + 55.6ei5^g^). 
The uniform solution is determined up to a free phase Lp which results from the orientation 
shift symmetry z ^^ ze^"^ of Eqn. ([2j). The positions of the pinwheels are fixed by the OD map 
and there are no translational degrees of freedom. In addition to these hexagonal pinwheel 
crystals (hPWCs) there exist also non-uniform solutions. Besides stripe-like solutions of z(x) 
with one dominant mode we find rhombic pinwheel crystals (rPWCs) Aj = Aj- = {A^ a, A) 
with a ^ A and distorted rhombic crystals Aj = (^1, ^2, ^3) , ^j- = (^3,^2,^1), both 
symmetric under rotation by 180 degree. We analytically calculated the stability properties 
of the uniform solution by linear stability analysis. The phase diagram for r^ <C r^, cf. 
Fig. [3](a) reveals a transition from rPWCs to hPWCs with increasing coupling strength e for 
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FIG. 3: (a) Phase diagram of the model Eqn. ([5|) for r^ ^ Tq. Vertical lines: stability range 
of OD hexagons, red solid line: stability border of hPWC, blue solid line: stability border of 
rPWC, blue dashed line: transition from rhombic to distorted rPWC. hPWC is ground state 
above the dashed red line, (b) hPWC (c) rPWC. OD contour lines: 10%, 50%, 90% 
contralateral eye dominance (black, gray, white), dashed line: unit cell. 

intermediate degrees of OD bias. For 7 < 7* or for 7 > 74 pinwheel-free orientation stripes 
are dynamically selected. For 7* < 7 < 74 and above a critical effective coupling strength 
ei5^ ^ 0.042, hPWCs are stable and become the energetic ground state above eB"^ ^ 0.12. 
Below eB^ ^ 0.065, rPWCs are stable leading to a bist ability region between rPWCs and 
hPWCs. We find in this region that rPWCs transform into distorted rPWCs above an ef- 
fective coupling strength of ei?^ ^ 0.033. Although rPWC are stable even in the uncoupled 
case they never become the energetic ground state. Thus for substantial bias towards one 
eye pinwheels are in fact stabilized and pinwheel rich solutions become ground states by 
inter-map coupling. 

The layouts of the main pinwheel rich solutions are shown in Fig. [3b, c. The hPWC contains 
6 pinwheels per unit cell and the pinwheel density i.e. the number of pinwheels per A^ jj] 
is given by p = 6 cos(7r/6) ^5.2. The rPWC has 4 pinwheels per unit cell and its pinwheel 
density is p = 4 cos(7r/6) ^ 3.5. One may expect that the energy term Eqn. ([3]) favors 



pinwheels to colocalize with OD extrema. For hPWCs three pinwheels of the same topolog- 
ical charge are in fact located at the extrema of the OD map. The other three however are 
located near OD borders. In case of the rhombic layout there is only one pinwheel at an OD 
extremum while the other three pinwheels are located at OD saddle-points which are also 
energetically favorable positions with respect to T. 

We tested whether these solutions and their stability ranges revealed for r^ <C r^ persist 
when the backreaction on the OD map is taken into account. To this end we solved the 
full field dynamics Eqn. ([2]) numerically using a fully implicit Krylov subspace algorithm 
with periodic boundary conditions on a 128x128 mesh with an aspect ratio of F = 22. In 
simulations we tracked the pinwheel density from t = to t = 10^ r~^, cf. Fig. [H In 
the uncoupled case (6 = 0), most of the patterns decay into a stripe solution and their 
pinwheel density drops towards zero. At small coupling strengths (e = 200) the pinwheel 
density converges either to zero (stripes), to values near 3.5 for the rPWC, or to approx- 
imate 5.2 for the hPWC. At high map coupling (e = 2000), pinwheel- free stripe patterns 
form neither from pinwheel rich nor from pinwheel free initial conditions. In this regime 
the dominant layout is the hPWC. However regions of hPWC layout can be inter-digitated 
with long lived rPWCs and stripe domains. Figure [4](d) shows the time course of the power 
P{t) = (|^(x, t)(i2/nP)x/(k(x, t)t/ip)x- The field Zth is obtained from solution of the ampli- 
tude equations Eqn. (|4j) while z^yn is the field obtained from the simulation. The amplitudes 
grow and saturate after t ^ r~^. When the amplitudes are saturated pattern selection 
starts. Quantitatively, we find that with backreaction the critical coupling strengths are 
slightly increased compared to their values in the limit r^ <C r^. 

Our analysis for the first time conclusively demonstrates that OD segregation can stabilize 
pinwheels, even if they are intrinsically unstable in the uncoupled dynamics of the OP map, 
raising the possibility that inter-map coupling is the mechanism of pinwheel stabilization in 
the visual cortex. Our results indicate that the overall dominance of one eye is important 
for the eflFectiveness of this mechanism. In this case, OD domains form a system of patches 
rather than stripes enabling the capture and stabilization of pinwheels by inter-map cou- 
pling. Studying a wide range of phenomenologically conceivable interaction energies we find 
that systems of OD stripes are in general not expected to stabilize pinwheel patterns. Inter- 
estingly, visual cortex around the time of early OP development is indeed dominated by one 
eye and has a pronounced patchy layout of OD domains supporting this notion [12]. Further 
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FIG. 4: Time evolution of the pinwheel density for r^ = 0.05, r^ = 0.25, 7 = 0.15. (a)-(c) 
Simulations in blue started from an identical set of 20 initial conditions. Red dashed lines: 
p = 4cos(7r/6) and p = 6cos(7r/6). (a)-(c) e = 0, 200, 2000. (c) OD and OP stripes as initial 
conditions (green), (d) Power of OP map, e = 0,200,2000 (red, blue, green). 

support comes from experiments in which the OD map was removed artificially, resulting in 
a significantly smoother OP map [10(]. Removal of the OD map, however, does apparently 
not completely destabilize pinwheels. This might refiect the infiuence of additional columnar 
systems of patchy layout like spatial frequency or direction columns that are expected to 
interact with the OP map in a similar fashion as OD columns. Interactions among multiple 
coupled maps may potentially also explain the non-crystalline spatial organization of OP 
maps in the visual cortex. The approach introduced here will be useful for further rigorous 
analyses of the interaction among multiple maps in cortical development. 
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